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Abstract. The surface-atom Casimir-Polder-Lifshitz force out of thermal 
equilibrium is investigated in the framework of macroscopic electrodynamics. 
Particular attention is devoted to its large distance limit that shows a new, 
stronger behaviour with respect to the equilibrium case. The frequency shift 
produced by the surface-atom force on the the center-of-mass oscillations of a 
harmonically trapped Bose-Einstein condensate and on the Bloch oscillations of 
an ultra-cold fermionic gas in an optical lattice are discussed for configurations 
out of thermal equilibrium. 



1. Introduction 

The electromagnetic force felt by a neutral atom near the surface of a substrate 
has been object of an intense investigation since the pioneering works by Casimir 
and Polder ^ and Lifshitz, Dzyaloshinskii and Pitaevskii [51 In addition to 
the fundamental character of the force, these studies U are presently motivated by 
the possibility of technological applications jSJ, by searching stronger constrains on 
hypothetical non Newtonian forces as well as its role in biological systems [7]. 

New perspectives to study such a force are opened by the recent development 
in storing and manipulating ultra-cold atoms. Indeed experimental and theoretical 
research has been recently focused on the forces acting on ultra-cold atomic gases due 
to the presence of a nearby surface. They include atomic beams 0, Bose-Einstein 
condensates [HI EH and degenerate Fermi gases ^21 ■ 

The surface-atom force at thermal equilibrium F'^^(T,z) can be in general 
separated in two parts 



The first one, Fq{z), is related to zero-point fluctuations (T = 0) of the electromagnetic 
field. At short distances z this force behaves like 1/z^ and is the analog of the van 
der Waals-London inter-atomic force. At larger distances the inclusion of relativistic 
retardation effects gives rise to the Casimir-Polder asymptotic behaviour |3] 



PACS numbers: 03.75.Kk, 67.40.Db, 77.22.-d, 78.20.-e 



F'^^iT,z)^Fo{z) + F:,^iT,z). 
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where ap and Eq are the static polarizabihty of the atom and the static dielectric 
function of the substrate respectively. The function (j){eo) ~ 1 is defined, for example, 
in [5]. The second contribution to the force, F^^{T,z), is due to the thermal 
fluctuations of the electromagnetic field. This contribution was first considered by 
Lifehitz who applied the theory of electromagnetic fluctuations developed by Rytov 
|14| . At large distances the thermal contribution approaches the so-called Lifshitz law 

W,z)_oo = -^^^^. (3) 

Such asymptotic behaviour is reached at distances larger than the thermal wavelength 
At = hc/ksT, corresponding to ~ 7.6 fim at room temperature. Thus it is the leading 
contribution to the total force. 

The Lifshitz force was originally evaluated at thermodynamic equilibrium. A 
non-trivial issue is the study of such a force out of thermal equilibrium, characterizing 
configurations where the temperature of the substrate Tg and environment Te, do 
not coincide. For instance in typical experiments with ultra-cold atomic gases the 
environment temperature is determined by the chamber containing the substrate and 
the trapped atoms. 

In this paper we describe the surface-atom force out of thermal equilibrium and 
how to recover its asymptotic behaviour at large distances. We assume that the 
radiation surrounding the atom is not able to populate its internal excited states 
which are assumed to be at energies hujat much higher than the thermal energy: 

ksTs, /cbTe < huat ■ (4) 

This condition is very well satisfied at ordinary temperatures (for example the first 
optical resonance of Rb atoms corresponds to 1.8 10''K). In the last part of the paper 
we analyze the eflfects of such a force on cold atoms, and in particular on the center- 
of-mass motion of a trapped Bose-Einstein condensate and on the Bloch oscillations 
of ultra-cold fermionic atoms in an optical lattice. 

2. Green-function formalism 

In the calculation of the surface-atom force the main ingredient is clearly the 
electromagnetic field and its sources. The latter, in our approach, are treated as point- 
like oscillating dipoles. Furthermore it is useful to write the fields using the Green's 
function formalism, the Green's function being the solution of the wave equation for 
a point-like source. Once this solution is known, the solution due to a general source 
can be obtained by the principle of linear superposition. The dyadic Green function 
G describing the electromagnetic field in surface optics (for isotropic, linear and non- 
magnetic media) is the solution of the equation 

V A V A G[lu; r, r'] - k^e{uj; r)G[w; r, r'] = ATTk'^I5{r - r'), (5) 

with the boundary conditions imposed by the geometry of the problem. In previous 
equation fc = w/c is the vacuum wavenumber, e(w;r) is the dielectric function and I 
is the identity dyad. Equation ij^l comes from the usual wave equation for the Fourier 
transformed electric field 

VA VAE[w;r] - fc2£(w;r) E[w;r] = 47rfc2p[w; r], (6) 

obtained from the macroscopic Maxwell equations in which the sources are described 
by the effective electric polarization field P[w;r] related to the electric current by 
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J[a;;r] = —iujP[uj;r]. The convolution of the solution obtained from Eq.© and the 
effective electric polarization gives the electric field at the observation point r 

E[u;;r]= jG[u;;v,r'] ■P[Lj;r'] dr'. (7) 

3. Surface-atom force 

Let us consider the atom described by its complex dielectric polarizability function 
a{u!) = a'{uj) + ia" {lo) in a vacuum half space Vi and placed at a distance z from the 
surface of the dielectric half space V2 ■ Let us choose an orthogonal coordinate system 
with the xy plane coincident with the interface and the z axis such that the dielectric 
occupies the region with z < and the vacuum the region with z > 0. The force 
acting on a neutral atom without a permanent electric dipolc moment is |15| 

F(r) = idr{t)VEr\v\t)\^ « iS^''{t)VEf{v\t)\^) + id^{t)VET\v',t)\^ , (8) 

where diS are the atomic electric dipole components, we have used the Einstein's 
summation convention for repeated indices and V' = Vr' . In Eq. the average is 
done with respect to the state of the atom and of the field and the lowest order in 
perturbation theory has been considered. The first term describes the (spontaneous 
and thermal) field fluctuations correlated with the induced dipole, and the second 
term involves (spontaneous and thermal) dipole fluctuations correlated to the field 
they induce. The induced electric dipole for the atom at the position r is 

d'"'i[cj] = a{w) E'°*[tj;r] « a{Lo) E"[tj;r], (9) 

where E^[a;;r] is the fluctuating field, and now a{ijj) is the atomic polarizability of 
the atom in an unbounded space. By modeling the atom as a point-like source dipole 
d{t) = d[aj]e~*"* at r, the corresponding polarization in the frequency domain is 
P[a;, r"] = d[uo] 5{y" — r), and the electric field at the position r' is 

E'"'i[cj;r'] = G[c^;r',r] •d*°*M « G[w;r',r] •d^'H. (10) 

Using Eq. and (|l()|l . the fluctuating dipole and field contributions to the surface- 
atom force © read 

drWV'i?^(r',t)[)= j I ^^e-(— ')*«(..) V'(i?f[c.;r]i?«t[,';r'])L, (11) 

df(t)V'4"<^(r',t)[)= J J ^^e-(— ')*V'G;,[c.;r',r][ (rff Hd«t[c.']) , (12) 
where the integrations are over the whole real frequency axis. 

4. Surface-atom force at thermal equilibrium 

At thermal equilibrium, in order to calculate the average values in 1)11(1 and 1(12(1 . 

it is possible to use the fluctuation dissipation theorem [141 [T?>j . One finds for the 
fluctuating dipoles 

d»Md;'M)^ ^-^^^^X!^ »"M. (13) 



and for the fluctuating flelds 

fltr„/..'i\ 



Ef[.-v] E^^W-A) ^ ImG.,[.;r,r']. (14) 
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After substituting the previous equalities into Ea. Hll|) and H12|) and using the 



reciprocity theorem Gy [a;;r, r'] 
equihbrium becomes 



Gji[Lu;r',r], the surface-atom force at thermal 



duj coth 



huj 
2fcBT 



Im 



a{Lu) dz Gii[uj] r, r'] 



(15) 



Because of the relation coth {nu/2kBT) = 1 + 2n {u/T), where n [oj/T) = (e'^'^/'^^T _ 
1)~^ is the Bose factor, one can separate in Ea. (|15|l the zero-point fluctuations 
contribution Fq{z) from the thermal contribution F°^{T,z). The latter term is the 
sum of two contributions arising from the two terms of Eq.©. The first one is due 
to the field fluctuations and it is linear in a' . The second one arises from the dipole 
fluctuations and it is linear in a" . As long as the condition Q is valid, the field 
fluctuations contribution is the leading term in F^^{T, z). 



5. Surface-atom force out of thermal equilibrium 

A first important investigation of the surface-atom force out of thermal equilibrium 
was carried out by Henkel et al. ^3] who calculated the force generated by a dielectric 
substrate at finite temperature by assuming that the environment temperature is zero. 
The principal motivation of that paper was the study of the force at short distances. 

In this section we analyze the general case of an atom placed in vacuum at a 
distance z from the flat surface of a substrate that we assume to be locally at thermal 
equilibrium at a temperature Tg which can be equal or different from the environment 
temperature Te, the global system being in or out of thermal equilibrium respectively, 
but in a stationary regime In this configuration it is relatively easy to describe 

the radiation produced by the flat substrate,while it is less trivial to describe the 
radiation coming from the environment. To face this problem we use the Lifshitz trick 
12] for which the vacuum half space is assumed to be a dielectric locally at thermal 
equilibrium with temperature Te, by introducing an infinitesimal imaginary part of 
its dielectric function. Using the fluctuation dissipation theorem and after integrating 
over an infinite volume the vacuum half space produces a radiation corresponding to 
the one that in a real systems is generated by the environment walls at Te. 

We refer to the substrate as to the half space 2 occupying the volume V2 with 
z < 0, with dielectric function £2{uj) — £2('^) + ^^2^'^) ^'^d in thermal equilibrium at 
the temperature Tg. The vacuum half space 1 instead occupies the volume Vi with 
z > and is characterized by a dielectric function £1(0;) = £'i{lo) -\- is'Kcu) and a 
temperature Te. Only after calculating the electric fields in this configuration we set 
ei{Lo) = 1. 

As well as for the thermal equilibrium case, the surface-atom force out of thermal 
equilibrium can be written as 

T-^(Ts, Te, z) - Foiz) + F^^l^s, Te, z), (16) 

where the thermal contribution F^^'^{Ts, Te, z), provided the condition Q is satisfied, 
is dominated by the thermal part of the fluctuating fields correlation (|ll|l only, as at 
thermal equilibrium! . 

The physical origin of the electromagnetic field is m] the random fluctuating 

J It is worth noticing that since zero-point fluctuations are not affected by condition jlj, in the 
calculation of the zero temperature force Fo{z) both dipole zero-point fluctuations 1131 and field 
zero-point fluctuations I14i are needed. 
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polarization field P[cj;r], whose correlations, at thermal equilibrium, are described 
by the fluctuation dissipation theorem 

/or jp*r ' S{uj ~uj')6(r -r')5ki he"{uj) 

Since the correlations of the source polarization field are local, the fluctuations of the 
sources at different points add up incoherently. Therefore we can assume that in the 
whole space the correlations of the sources are given by equation (|17() , valid for source 
dipoles in the the two half-spaces assumed to be locally at thermal equilibrium at two 
different temperatures [211. In order to calculate the field correlation function 
we express the electromagnetic field in terms of its source polarization field via Eq. 
O, and using the Eq. (|17|l we write the thermal part of the surface-atom force out 
of thermal equilibrium as 

(Ts, Te, z) = (Ts, 0, z) + F^iO, ^e, z), (18) 

where the first thermal contribution 

4'(^)Re \a{cj) G.fc[u;; r, r']d.G*,[u;; r, r'] d^r'" 

^^tr(Ts,o,z) = ^y^ dco ^ e^w^BTs _ 1 (19) 

arises from the sources in the substrate V2, while the second one 

^^tr(0,TE,z) = ^y^ do. ^ ,n./.eT. _ 1 (20) 

is produced by the sources in the vacuum half space Vi §. It is possible to show that 
the sum of Ea. lfT^ and (jJOJ, at the same temperature, reproduce the thermal part of 
the force at thermal equilibrium JS| ■ Indeed it is possible to apply to the whole space 
the Green's functions property (see, for example, |19|) 



j dr e"(r,Lj) Gifc[w; ri, r]G*fe [w; r2, r] = 47rlm Gy [cj; ri, r2]. 



(21) 



where the integration is on the volume such that on its surface the Green function is 
zero. Than we can express the complete surface-atom force out of thermal equilibrium 
in the convenient form 

F-'i(Ts, Te, z) = F'^^T^, z) + F^^\Ts,0, z) - F^^\T^, 0, z), (22) 

where the equilibrium force F'^'^{T,z) is given by ^ and F^^'^{T,0, z) is defined by 
Ea. H19l) . Consistently with the assumption Q), in deriving the thermal part of Ea. l22|l 
we ignored terms proportional to the imaginary part of the atomic polarizability. For 
the same reason the wind contribution in Ea. H19|l and H20() . related to a", can be 
ignored and the real part a'(w), corresponding to the dispersive contribution, can be 
replaced with its static {lo = 0) value ao ■ In this non-absorbing condition the force of 
Ea. H18|) can be also written in the form F^^'^{Ts,Te, z) = 47rao9zC^Ei('7E, Ts, z) where 
Uea = {E{z,t)^') /Btt is the thermal component of the electric energy density at the 

§ The Green function Gik then reduces respectively to its transmitted component in Eg. 1191 1211 
and to its incident, reflected and local component in Ea. t2UI 1181 . 
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z [ilm] 

Figure 1. Surface-atom force F^^'^'^(z) calculated from Ea. l22i . for different 
thermal configurations. 



atom position. 

After some lengthy algebra we find for Eq. (|19|) the relevant result 



poo 4 poo 



xV|eH-92| + (e'H-g2) 



^£(w) - + ^1 - g2 



r + 



^/s{uj) - q^ + £{uj)yjl - q' 



(23) 



where we introduced the dimensionless variable q — Kc/lj^ with K the modulus of the 
electromagnetic wave- vector component parallel to the interface, and s{ijj) = S2{'^)- 
In Figure ^ we show the explicit results for the full force obtained from Eq. H22|l 
as a function of the distance from the surface for different choices of Ts and Te. 
The calculations have been performed for a sapphire substrate (eq — 9.41) and for 
*^Rb atoms (ao = 47.3 10"^'* cm^). For F'^'^{T,z) we have used the predictions of 
P]. The figure clearly shows that the thermal effects out of equilibrium are sizable 
(solid lines), thereby providing promising perspectives for future measurements of 
the surface-atom force at large distances. To increase the attractive nature of the 
force it is much more convenient to heat the substrate by keeping the environment 
at room temperature (lower solid line) rather than heating the whole system (dashed 
line). When Tg < Tg (upper solid line) the force exhibits a characteristic change of 
sign reflecting a repulsive nature at large distances (see also discussion below). At 
short distances the thermal correction to the force becomes smaller and smaller and is 
determined by the temperature of the substrate. The new effects are visible already 
at distances z = 4-^ 7/im, where experiments are now becoming available |10| . 
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6. New asymptotic large distance limit 

In this section we discuss in details the large z behaviour )17| of the out of equilibrium 
force (|221l- After the substitution q'^ - 1 ^ t"^, Eq.|(23|l becomes 

KriTA^)--'-^ /f du. ^,^/;,_ J °°dtt^e-^-/-m.), (24) 
where 

/ 

f{t,iu) = ^\e{u;)-l-t^\ + is'iu)-l-t^) 



y/e{Lu) -l-t^ +it 



{2t'^ + l) (1 + t^ + \e{uj) - 1 - t 



21 



yeH^T^^ + is{Lu)t 



(25) 



F;r(T,0,z)_eo^-^ / du.—j£—f{^). (26) 



Due to the presence of the exponential e^^^*'^/'^ in Ea. H24|l . it is possible to show that 
only the region i <C 1 contribute to the large z behaviour of the force that in such 
limit exhibits the non trivial asymptotic behaviour 

This force exhibits a slower 1/z^ decay with respect to the one holding at thermal 
equilibrium where it decays like l/z** (see Eq. ©). In the above equation we have 
introduced the low t expansion of Eq. 125() 

/ (c) = y^leico) - 1| + [e'iuj) 1] ^ l^l^)"!;, ^' ■ (27) 

Result H26|l and (|27|l provide the large distance behaviour {z — > oo) of the force (|23|l 
where the only assumption made was the condition Q . Due to the presence of the Bose 
factor the force 126|l depends on the optical properties of the substrate at frequencies 
of the order of ~ k^T/fi. 

For temperatures much smaller than HuJc/kB, where coc is the lowest characteristic 
frequency of the dielectric substrate, only the static value of the dielectric function is 
relevant and so we can replace f{t,Lu) with its low frequency limit in Ea. H24|) . In this 
limit f{t,uj 0) is different from zero only for < t < -^eo — 1, and after the t 1 
expansion Ea. l|24(l becomes 



where we performed the change of variable u = 2ztuj/c and replaced f{t,uj — > 0) with 
its i ^ 1 expansion \/2(eo + l)/\/^o — 1- For 

where At — fiLj/k-B,T is the thermal wavelength, we can extend the upper limit 
of integration on u to +oo and so we obtain that the force (|22|l felt by the atom 
approaches the asymptotic behaviour 

F (Ts,TE,z).^oo--g -3^^ (30) 
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Result (|30ll holds at low temperature with respect to the first dielectric function 
resonance (T ^ TiWc/^b) and at distances satisfying the condition (|29|l calculated 
at the relevant temperatures Tg and Te- Ea. (|30|l shows that, at large distances, the 
new force is attractive or repulsive depending on whether the substrate temperature 
is higher or smaller than the environment one. Furthermore, it exhibits a stronger 
temperature dependence with respect to equilibrium force ((SJ, contains explicitly the 
Planck constant and has a l/z"^ distance dependence! | . 

The new dependence of F^°'^{T, 0, z) on temperature and distance can be 
physically understood by noticing that the main contribution to the z dependent 
part of the electric energy Uei arises from t <C 1. Such values of t correspond to the 
component of the black-body radiation impinging on the surface from the dielectric 
side in a small interval of angles, of order of (At /z)'^, near the angle of total reflection. 
This radiation creates slowly damping evanescent waves in vacuum. As a result 
pneqf^rp^ 0, z) tums out to be, in accordance with Ea. l|Sn|) . of order of —{aQX'^/z^)UBB, 
where Ubb oc T'* is the energy density of the black-body radiation. 

Equation (|30|) holds for a dielectric substrate where Eq is finite. If we want to 
find the large distance limit for a metal we should use equation H26() . In the limit 
of small values of T we can use the Drude model. As only frequencies us ~ fceT/Ti 
contribute, one can substitute in Eq.(j2Zjl = 47rcr/a; 3> 1, the real part e'{uj) 

remaining finite as w — + 0. Than one finds / (cj) y^e"{uj) = 2-^/ na/uj, where a is 
the electric conductivity, so that for a Drude metal Ea. H30() is replaced by 

F (Ts,Te,z)__ = , (31) 



where C(3/2) ^ 2.61 is the Riemann function. It is easy to show that Ea.l|31|l is valid 
at the condition 

z > Tc■/^c^/A^/{kBT f'^■ (32) 



7. Effects of the surface-atom force on ultra-cold atoms 

Ultra-cold gases can provide a useful probe of the surface-atom force. A mechanical 
tool sensitive to the gradient of the surface-atom force is in fact the frequency shift 
of the center-of-mass oscillation of a trapped Bose-Einstein condensate EI • On 
the other hand, experiments based on Bloch oscillations are interferometric tools 
sensitive to the force itself Ej- Finally one could also think to interference 
experiments involving the macroscopic phase of Bose-Einstein condensates in a double 
well potential |^ EH]. The position of the corresponding interference fringes are 
sensitive to the surface-atom potential. In the last part of this paper we discuss the 
first two above mentioned experiments. 



7.1. Effects on the collective oscillations of a trapped BEC 

Bose-Einstein condensed gases I^HI are very dilute, ultra-cold samples characterized 
by unique properties of coherence and superfluidity. The study of the collective 
oscillations '571 of a Bose-Einstein condensate provides a useful probe of the surface- 
atom potential. In fact it is possible to measure with great accuracy the frequency 

I Instead of calculating the asymptotic behaviour 1301 of the force from the general equation 121^1 . 
it is possible to produce a more direct derivation assuming from the very beginning that one can 
neglect absorption and dispersion of the dielectric function of the substrate 1221 . 
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Figure 2. (a) Relative frequency shift I33i of the center-of-mass oscillation of a 
BE condensate {Rz = 2.5/im, liJz/2it = 220Hz) and (b) relative shift of the Bloch 
oscillation period of a degenerate Fermi gas 1341 , out of thermal equilibrium I16i . 



of the center-of-mass motion lucm of a condensate. For a harmonically trapped 
condensate the frequency ujcm corresponds to the harmonic trap frequency Uz , where 
z is the direction of the oscillations. Thus if a BEC in a harmonic trap is placed at 
distance d from the surface of a substrate, the surface-atom potential Vs-a{z) perturbs 
the trap potential and produces a shift in locm- In the limit of small oscillations (for 
a complete analysis see -Q ), such a shift is 

2 r+R, 

^CM^^l + - / dzn'oiz)dys-.{z + d), (33) 

where UqIz) is ID column density of the gas (density integrated over the directions 
perpendicular to the direction of oscillation) and Rz is the Thomas-Fermi radius in 
the z direction^. Therefore measuring wcm it is possible to extract the surface-atom 
potential V;_a(^) PITU]. 

In fig. Ola) we plotted, as a function of the surface-condensate separation d, the 
relative frequency shift Aujqm/^z = {^z —^cm)/'^z for the center-of-mass oscillations 
of a *^Rb condensate close to a sapphire substrate. In such a calculation we used the 
surface-atom potential corresponding to the force 1)22(1 . 

7.2. Effects on Bloch oscillations in Fermi gases 

When an external force Fcxt is applied to a particle trapped in a periodical potential, 
the particle undergoes oscillations in momentum space (the Bloch oscillations) . During 
this oscillations the particle quasi- momentum q evolves according to hq = i^cxt- This 
is what happens for example in a sample of ultra-cold atoms trapped in a ID optical 
lattice aligned along the vertical direction. Bloch oscillations produced by the effect 
of the gravity force Fq = mg have a period Tb — AttU /mg\ where A is the lattice 
wave-length and g is gravity acceleration. If now a surface is brought close to the 
atomic sample the additional surface-atom force Fs-a.{z) perturbs the gravitational 
potential and affects the dynamics of the Bloch oscillations (for a complete analysis 



^ For a Bose-Einstein condensate in its ground state the ID column density is easily evaluated in 
the so called Thomas-Fermi approximation where nQ(z) = 15(1 — / lij^)^ / IQRz I2bl . 
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see ^1 ESI)- In particular it produces a shift of the period Tb. In figure [2Ib) we 
plotted the relative shift ATb/Tb for different thermal configuration as the distance 
d between the center of a cloud of fermionic atoms (ao = 4.3 10~^'^cm^) and the 
surface of a sapphire substrate is varied. We used also the approximation of a small 
cloud of Fermi atoms, for which 

^ = (34) 

Tb mg 

In the range of distances plotted in fig. EJb) this approximation provides results in 
good agreement with the exact calculation [l2j that takes into account real experi- 
mental parameters of the gas. 

It is worth noticing that both effects of the surface-atom force out of thermal 
equilibrium described in the last section, and plotted in Figures 121 a) and|2Ib), are in 
the domain of the present experimental accuracy JOl HHj ■ 
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